Magnetic properties of a strongly correlated system on the Bethe lattice 
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We study the influence of an external magnetic field h on the phase diagram of a system of Fermi 
particles living on the sites of a Bethe lattice with coordination number z and interacting through on- 
site U and nearest-neighbor V interactions. This is a physical realization of the extended Hubbard 
model in the narrow-band limit. Our results establish that the magnetic field may dramatically affect 
the critical temperature below which a long-range charge ordered phase is observed, as well as the 
behavior of physical quantities, inducing, for instance, magnetization plateaus in the magnetization 
curves. Relevant thermodynamic quantities - such as the specific heat and the susceptibility - are 
also investigated at finite temperature by varying the on-site potential, the particle density and the 
magnetic field. 
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I 

C 

o 
o 



> 

o 
o 



X 



Statistical models on the Bethe lattice are of consid- 
erable interest since they admit a direct analytical ap- 
proach for a number of problems that may be other- 
wise intractable on Euclidean lattices. Due to the pe- 
culiar structure of the lattice, several interesting physi- 
cal problems involving interactions are exactly solvable 
when defined on the Bethe lattice There is a gen- 
eral interest in the study of models defined on such lat- 
tices which goes beyond physics. Being a dendrimer of 
infinite generation, the Bethe lattice is attractive both 
for basic and applied interdisciplinary research involving 
chemistry, physics, biology, pharmaceutics, and medicine 
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In a recent paper, we provided a comprehensive and 
systematic analysis of the extended Hubbard model in 
the narrow-band (atomic) limit (AL-EHM) on the Bethe 
lattice with arbitrary coordination number z [3j . Within 
the Green's function and equations of motion formalism, 
we exactly solved the model and, by considering relevant 
physical quantities in the whole space of the model pa- 
rameters, we investigated the finite temperature phase 
diagram, for both attractive and repulsive on-site and 
intersite interactions. The aim of the present paper is 
twofold. First, we would like to further develop our pre- 
vious work, by extending it to a more general situation 
in which a finite magnetization may be induced by an 
external magnetic field. Secondly, the AL-EHM on the 
Bethe lattice exhibits interesting features at low temper- 
atures such as the existence of magnetic plateaus. Here, 
we address the problem of determining the effect on the 
phase diagram and on the behavior of several thermody- 
namic quantities of the presence of a uniform magnetic 
field, introduced through a Zeeman term. We study the 
properties of the system as functions of the external pa- 
rameters n, T/V, U/V and h/V, allowing for the on-site 
interaction U to be both repulsive and attractive. In 
fact, the parameter U can represent the effective interac- 
tion coupling taking into account also other interactions. 
Throughout the paper, we consider a repulsive intersite 
interaction V and we set V = 1 as the unit of energy, tak- 



ing the Boltzmann's constant fc# = 1. In the absence of a 
magnetic field, the phase diagram in the plane (n, T) ex- 
hibits a transition line along the temperature axis, below 
which translational invariance is broken Q. The Bethe 
lattice effectively splits in two sublattices with different 
thermodynamic properties. As a result, a charge ordered 
(CO) phase, characterized by a different distribution of 
the electrons in alternating shells, is established for perti- 
nent values of the particle density. A finite magnetic field 
can have dramatic effect on the phase diagram: for in- 
stancejwhen U < 0, the reentrant behavior observed for 
h = 3] disappears and one finds, above a certain value 
of h, a CO phase for n > 1/z. This should be compared 
with the case h = 0, where the CO phase is observed 
only for n > 2/z when U < 3] . The characteristic lobe 
structure, exhibited by the critical temperature as a func- 
tion of the particle density, shrinks by augmenting the 
magnetic field. By further increasing h, the CO phase is 
suppressed in the particle density range 0.75 < n < 0.85 
and 1.15 < n < 1.25 at h = h T = zV/2 - U/2; the lobe 
splits in three separated lobes centered around n = 0.5, 
1, and 1.5, respectively. The central lobe eventually van- 
ishes by increasing h. A similar behavior is observed also 
for U > 0. 

The magnetic properties of the system depend on the 
value of the particle density and of the on-site poten- 
tial. At low temperatures, and for attractive on-site in- 
teractions, the magnetic field does not play any role if 
its intensity is h < \U\/2: the ground state is a col- 
lection of shells with doubly occupied sites (doublons) 
surrounded by empty shells. The magnetic energy is not 
strong enough to break the doublons. On the other hand, 
for strong repulsive on-site interactions, it is sufficient a 
small nonzero value of the magnetic field to have a finite 
magnetization. In the intermediate region, the competi- 
tion among U, V and h determines the phase structure. 
For all values of the particle density, one finds a critical 
value of the magnetic field h s - dramatically depending 
on U and n - above which the ground state is param- 
agnetic. In this state, every occupied site contains one 
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and only one electron, aligned along the direction of h. 
For strong repulsive on-site interactions, h s = 0, i.e., the 
spin are polarized as soon as the magnetic field is turned 
on. Furthermore, for attractive on-site interactions and 
0.5 < n < 1, one observes the existence of two criti- 
cal fields, namely: h c , up to which no magnetization is 
observed, and h s , marking the beginning of full polariza- 
tion. This is analogous to the finite field behavior of the 
S = 1 Haldane chain [J] . 

The addition of a homogeneous magnetic field does not 
dramatically modify the framework of calculation given 
in Ref. [H, provided one takes into account the break- 
down of the spin rotational invariance. For the sake of 
comprehensiveness, in the next section, we briefly report 
the analysis leading to the exact solution of the AL-EHM 
on the Bethe lattice in the presence of a magnetic field. 
Then, we investigate the phase diagram in the space n, 
T/V, U/V by varying h and we observe that, above a crit- 
ical value of h, the CO region shrinks due to the presence 
of a finite magnetic field. We also analyze the magnetic 
properties of the system and find magnetic plateaus in 
the magnetization curve. Finally, last section is devoted 
to our conclusions. 



I. EXACTLY SOLVABLE MODEL 

When defined on the Bethe lattice with coordination 
number z, the narrow-band limit of the extended Hub- 
bard model, in the presence of an external homogeneous 
magnetic field, can be described by the following Hamil- 
tonian: 

z 

H = -fin(Q) + UD(Q) - hn 3 (Q) + ^H^ , (1) 

P =i 

is the Hamiltonian of the p-th sub-tree rooted at 
the central site (0) and can be written as 



double occupancy operators at site i, respectively, n 3 (i) 
is the third component of the spin density operator, also 
called the electronic Zeeman term, 



= nf(i) - n^i) = c\(i)c^(i) - c[(i) Ci (i). (4) 

Here we do not consider the orbital interaction with 
the magnetic field. As usual, n a (i) — c\(i)c a {i) with 
a = {t)41) where c a (i) (c\(i)) is the fermionic annihila- 
tion (creation) operator of an electron of spin a at site i, 
satisfying canonical anticommutation relations. We use 
the Heisenberg picture: i = (i,i), where i stands for the 
lattice vector R4. 

The exact solution of the model can be obtained by 
using the equations of motion approach in the context 
of the composite operator method || , which is based on 
the choice of a convenient operatorial basis. For our pur- 
poses, the suitable field operators are the Hubbard oper- 
ators, = [1 — n{i)]c a {i) and rj a {i) = n(z)c CT (i), which 
satisfy the equations of motion: 
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d 

l Jf[. T ]<?( t ) = "(A* - U + <7h)r] a (i) + zVr] a (i)n a (i). 



(5) 



Hereafter, for a generic operator $ (?) we shall use the 
notation $"(z) = J2p=i ®(hP)/ z i (hP) being the first 
nearest-neighbors of the site i. The Heisenberg equa- 
tions (O contain the higher-order nonlocal operators 
5cr(i)n Q (i) and r] <T (i)n a (i). By taking time derivatives 
of the latter, higher-order operators are generated. This 
process may be continued and an infinite hierarchy of 
field operators is created. However, since the number 
n(i) and the double occupancy D(i) operators satisfy the 
following algebra 



n p (i) — n(i) + a p D(i), 
D p (i) = £>(i), 
n p (i)D(i) = 2D(i)+a p D(i), 



(6) 



where p > 1 and a„ = 2 P — 2, it is straightforward to 



H {v) = ~nn(p)+UD{p)-fni 3 (p)+Vn(0)n(p)+Y^ ff (p ' m) . establish the following recursion rule 



m— 1 



(2) 

Here, (p) (p = 1, . . . z) are the nearest- neighbor sites of 
(0), also termed the first shell. ij(P< m ) describes the m-th 
sub-tree rooted at the site (p): 



2z 



(7) 



jj{p,m) = _^ n (p )m ) + U D(p,m) - hn 3 (p,m) 

z-l 

+ Vn(p) n(p, to) + H^" 1 '^ . 

9=1 



(3) 



(p, m) (to = 1, ... z — 1) and (0) are the nearest-neighbors 
of the site (p). H^ p ' m ^ is the Hamiltonian of the q- 
th sub-tree rooted at the site (p,m). The process may 
be continued indefinitely. U and V are the strengths 
of the local and intersite interactions, respectively; /1 is 
the chemical potential, n(i) = n^(i) + n^(i) and D{i) = 
n^(i)ni(i) = n (i) [n (i) — 1] /2 are the charge density and 



which allows one to write the higher-power expressions 
of the operator n a (i) in terms of the first 2z powers. 

The coefficients iJJ' are rational numbers, satisfying the 

relations J2m=i ^ = 1 and ^™ = ^m,fc (k = l,...,2z) 
0, @] ■ The recursion relation ([7]) allows one to close the 
hierarchy of equations of motion. As a result, a complete 
set of eigenoperators of the Hamiltonian ([1]) can be found. 
To this end, one defines the composite field operator 



(8) 
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where 



\t*(i)[n a {i)] 2 'J 

( nS) \ 



(9) 



With respect to the case of zero magnetic field @, the 
degrees of freedom have doubled, since one has taken into 
account the two nonequivalent directions of the spin. By 
exploiting the algebraic properties of the operators n(i) 
and D(i), and the recursion rule (O, it is easy to show 



that the fields ip^ (») and are eigenoperators of 

the Hamiltonian Q @: 



e^' and e' 17 - 1 are the energy matrices: 
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(10) 



, (11) 



where e& and ei^ are the (2z + 1) x (2z + 1) matrices 
given by 
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U-n-ah 
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zV 



(13) 



U - fi - ah + zVA 



2z+l 
2z I 



The eigenvalues of the matrices e£P and are 



= -/i- ff ft + (p-i)v; 

E($ = U-»-ah+(p- 1)V, 



(14) 



with p=l,...,2z + l. The Hamiltonian has now been 
formally solved since, for any coordination number of the 
underlying Bethe lattice, one has found a closed set of 
eigenoperators and eigenenergies. As a result, one may 
solve the model and compute observable quantities. 

Since the addition of a homogeneous magnetic field 
does not dramatically modify the framework of calcula- 
tion given in Ref. Q , here we report only some details of 
the calculations and refer the interested reader to Ref. @ 
for a comprehensive analysis. Upon splitting the Hamil- 



tonian (JT|) as 



H = H, 



H 



(i) 



H { j i] = zVn(i)n a (i), 



(15) 



it is immediate to notice that, with respect to the case 

(i) 

h = 0, only is modified by the presence of the mag- 
netic field h. Therefore, all the calculations relative to 
Hj 1 ^ are not modified; the changes induced by the pres- 
ence of h will concern only the calculations involving Ho. 
In particular, the statistical average of any operator O 
can be expressed as 



(O) 



(Oe-^' ) ) , i 



1 )0,i 



(16) 



The symbol (• • • )o.i stands for the thermal average with 
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respect to the reduced Hamiltonian H^: i.e., (• • • )o,i — 

Tr{- ■ ■ e- pH ^ }/Tr{e- f3H o ) }. Equation {Tj| allows us to 
express the thermal averages with respect to the complete 
Hamiltonian H in terms of thermal averages with respect 
to the reduced Hamiltonian Ho, which describes a system 
where the original lattice has been reduced to the site i 
and to z unconnected sublattices. As a consequence, in 
the iJo-representation, correlation functions connecting 
sites belonging to disconnected sublattices can be decou- 
pled. Within this scheme, all the local correlators neces- 
sary to compute the Green's functions [3] can be written 
as functions of the parameters 



1 z 



P =i 



J5 = ^(i»W = {B%))o,i, 



(17) 



p=i 



in terms of which one may find a solution of the model. 
In the above equation, i p (p = 1, . . . , z) is an arbitrary 
neighboring site of i. 

A repulsive intersite interaction disfavors the occupa- 
tion of neighboring sites. At low temperatures, this may 
lead to a CO phase characterized by a nonhomogeneous 
distribution of the electrons in alternating shells Q . In 
order to capture this phase, we shall divide the lattice 
into two sublattices: A contains the central point (0) and 
the even shells, the sublattice B contains the odd shells. 
Then, one requires the following boundary condition to 
hold: 



("(<)> 



n A 
n B 



i G A, 
i e B. 



(18a) 



and 



{ip)) = - [kXi + f(l + k 2 ) K {X, + 2aY l ) F* 



+ kgK 2 [X l + 2dY i )Gr 1 ] , 



kY, 



(D(i p )) = -f [1 + /(l + k^K'Fr 1 + gkK^Gr 1 } 



(21) 



where / = e^ M , g = e' 3 ^ 2 ^ u \ and k = e l3h , and we used 
the definitions 



Fi = l + aXi + a z Yi, 



G i = l + dX l + d 2 Y l , (22) 
F=k + f(l + k 2 )F: + kgGt, 

K - 1, and d = K 2 - 1. Simi- 
larly, also the third component of the spin density can be 
written as a function of the parameters Xi and Yi : 



(n 3 (i)) 



f(l-k 2 )F? 



k + /(l + fc 2 )^ +gkG z t 



It is not difficult to show that the magnetization m — 
(713(1)) can be written in terms of the particle density and 
double occupancy as: (713(2)) = tanh(/3/i)[n(i) — 2D(i)]. 
Equations (fl9)) . together with Eq. (|1 8b[) - which fixes 
the chemical potential /i - constitute a system of coupled 
equations allowing us to ascertain the five parameters fj,, 
Xai Xb, Ya, and Yb in terms of the external parameters 
of the model, namely: n, h, U, V, and T. Once these 
quantities are known, all the properties of the model can 
be computed. 



II. PHASE DIAGRAM AND MAGNETIC 

-(n A +n B )- (18b) PROPERTIES 



Let us take two distinct sites i € A and j £ B. We 
require that the expectation values of the particle den- 
sity and of the double occupancy operators at the site i 
are equal to the ones of the neighboring sites of j and 
viceversa. As a consequence, the number of unknown 



correlators is four: the parameters Xi, Xj, 
determined by the equations 



Yi, Yi are 



(n(i)> = (n(j p )) 
(D®) = (D(jp)) 



(n(i p )) = (n(j)), 
(D(i p )) = (D(J)). 



(19) 



After lengthy but straightforward calculations, one can 
find analytical expressions for the average values of the 
particle density and of the double occupation operators 
in terms of the parameters Xi and Yi, namely: 



(n(i)> = 
(D(i)) 



f(l + k 2 )Ff +2gkGf 
k + f(l + k 2 )F? +gkG\ 

gkG? 

k + f(l + k 2 )F?+gkGl 



(20) 



In this section, we derive the phase diagram by numer- 
ically solving the set of equations (IT91 . In the absence of 
a magnetic field, we find regions of the (U, n, T) 3D space 
characterized by a spontaneous breakdown of transla- 
tional invariance Q. In these regions, the population of 
the two sublattices A and B is not equivalent: the system 
has entered a finite temperature long-range CO phase. 
Upon decreasing the temperature, the distribution of the 
electrons becomes more inhomogeneous. In the presence 
of a magnetic field, the phase structure is determined by 
the three competing terms of the Hamiltonian: the re- 
pulsive intersite potential (disfavoring the occupation of 
neighboring sites), the magnetic field (aligning the spins 
along its direction, disfavoring thus double occupancy) 
and the on-site potential, which can be either attractive 
or repulsive. The competition among these terms may 
affect the transition temperature. To understand the ef- 
fect of a magnetic field on the critical region, in Figs. Q] 
we plot the phase diagram at constant U for several value 
of ft.. In Figs. [T]we consider the full range of variation of 
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0.25 0.50 0.75 1.00 1.25 1.50 1.75 




FIG. 1: Phase diagram in the plane (T,n) for V = 1 and 
z — 3 and several values of h, and for (a) U = — 1; (b) 17 = 1. 



the particle density < n < 2, although it would be suf- 
ficient, owing to the particle- hole symmetry, to explore 
just the interval [0, 1]. 

In the plane (n,T), for attractive U, the CO phase 
is observed in a interval An which varies with the mag- 
netic field. As it is shown in Fig. [T^,, at h = and 
in the limit T — > 0, a complete CO state is established 
in the region 2/z < n < 2(z — l)/z. An first increases 
with T, then decreases vanishing at n = 1, where the 
maximum critical temperature is reached; a reentrant 
behavior characterizes this region. Upon turning on a 
finite magnetic field, An increases by increasing h and 
the reentrant behavior is lost. For h > \U\, a CO state 
is established in the region 1/z < n < (2z — l)/z (in the 
limit T — > 0). The phase diagram still presents a single 
lobe structure centered at n = 1, although the height of 
the lobe has decreased. There exists a critical value of 
the magnetic field Kt — zV/2 — U/2 above which one 
observes the formation of three separated lobes centered 
around n = 0.5, n = 1 and n = f.5, respectively. The 
transition to the CO phase is suppressed in the range 
0.75 < n < 0.85 and 1. 15 < n < f.25 (for z = 3). By 
further increasing h, the central lobe shrinks and even- 
tually disappears. As one can infer from Fig. [Ifc>, the 
picture is very similar for repulsive on-site interaction. 
The phase diagram presents a single lobe structure cen- 
tered at n = 1: a CO phase is observed below the critical 
temperature in the range 1/z < n < (2z — I)/z, which 
does not depend on h. At the same value of the magnetic 
field fir = zV/2 — U/2, the CO phase is observed inside 
the three separated lobes as for the case U < 0. When 
plotted as a function of the magnetic field, the critical 
temperature shows a decreasing behavior, which is more 
pronounced in the neighborhood of n = I . Moreover, as 
one can also notice from Figs. [TJ according to the value 
of the particle density, two situation can occur: either 
the transition temperature is finite for all values of h or 
it vanishes at same critical value. In Fig. [2^,, we plot the 
transition temperature T c as a function of the magnetic 
field at n = 0.75, z — 3 and for different values of the 
on-site potential. The transition temperature decreases 
by increasing U and vanishes at Kt = zV/2 — U/2. For 
U = 3V, fir = 0: a strong on-site repulsion inhibits the 
CO phase, as already noticed in Ref. @. The study of 
the specific heat also enlightens the influence of the mag- 




FIG. 2: (a) The critical temperature T c as a function of the 
magnetic field h at n = 0.75 for V = 1 and z = 3, and 
several values of U. (b) The specific heat as a function of the 
temperature at n = 0.75, for V = 1, U = — f , and z = 3, and 
several values of h. 



netic field on the thermodynamic behavior of the system. 
The specific heat is given by C = dE/dT, where the in- 
ternal energy E can be computed as the thermal average 
of the Hamiltonian (p} . As an example of the character- 
istic behavior of the specific heat by varying h, in Fig. [2Jd 
we plot C as a function of the temperature at U = — 1, 
2 = 3 and n = 0.75, for several values of the magnetic 
field. For h — and at low temperatures, the system 
is in a CO phase and exhibits a phase transition at T c 
to a homogeneous phase. The specific heat exhibits a 
peak at T\ = T c - due to the phase transition - and an- 
other peak Ti at low temperatures which vanishes as h 
increases. For finite magnetic fields, the position of the 
peak Ti and the relative height decrease by augmenting 
h. For h — hr, the transition is suppressed and, corre- 
spondingly, the peak T\ disappears. 

The competition among the magnetic field and the on- 
site and intersite potentials gives rise to the formation 
of plateaus in the magnetization curves. By increasing 
the magnetic field, one observes plateaus whose starting 
points depend on the particle density, as well as on the 
on-site potential: one identifies two critical values of the 
magnetic field. The nonzero magnetization can either 
begin from h — or from a finite field. h c denotes the 
starting point of a nonzero magnetization, whereas h s 
denotes the value of the magnetic field when it reaches 
saturation. The results for the magnetization m(h) are 
shown in Fig. [31 where U c is the critical value of the on- 
site potential separating the different observed behaviors; 
for z — 3, one finds U c f» 2.8. Making reference to Fig. [3J 
for the different regions of n and U , one has: h c = h s = 
\U\/2 (Fig. Oi), h c = \U\/2 and h s = zV/2+\U\/2 (Fig. 
Eb), h c = and h s = zV/2 - U/2 (Fig. [3b), and h c = 
h s = (Fig. [3b!) . The so-called metamagnetic behavior 
is clearly seen: at low temperatures the magnetization 
begins to show a typical S'-shape which becomes more 
pronounced by further decreasing the temperature. At 
T = one, two or three plateaus (m = 0, m = 1 — n 
and m — n) are observed, according to the values of 
the external parameters. These results are similar to the 
ones obtained in one- dimensional AL-EHM and spin-1 
antiferromagnetic Ising chain with single-ion anisotropy 



6 




■ T=0.2 

- 7=0. 1 

- 7=0.05 



o<u<u c 

0.5<n<1 



■7=0.2 
■7=0. 1 

- 7=0.05 



^ h 




(<& 




m = n 








7=0 2 

7=0. 1 

7=0.05 

7=0 




j! ■■' v=i 

V ' z=3 

V, u>o 

;,'■■ 0<n<0.5 


V=J 
z=3 
U> U 
0.5<n< 1 



o h , h h 

FIG. 3: The magnetization m as a function of the magnetic 
field h at T = 0. 



EH El- 

In Figs. H^-b we plot the spin susceptibility as a func- 
tion of the magnetic field at T = 0.1 for several val- 
ues of U (both attractive and repulsive), z — 3 and for 
n = 0.75. In the limit T —> 0, the spin susceptibility 
diverges in correspondence of the values h cr n at which 
the system moves from one magnetization plateau to the 
other. For low values of the magnetic field and attrac- 
tive on-site interactions - corresponding to Fig. [4ji - the 
spin susceptibility vanishes at low temperatures for all 
values of the filling: in the limit T — > all electrons are 
paired and no alignment of the spin is possible. By in- 
creasing h, the magnetic excitations break some of the 
doublons inducing a finite magnetization: Xs has a peak, 
then decreases, the system having entered the successive 
magnetic plateau. If 0.5 < n < 1 then another peak is 
observed, corresponding to the second jump of the mag- 
netization when h reaches the saturated value h s . On 
the other hand, for repulsive on-site interactions, a very 
small magnetic field induces a finite magnetization (with 
the exception of n = 1 when U < U c ): Xs has a maxi- 
mum at h = and then decreases by augmenting h, un- 
less another transition line is encountered, as it happens 
for 0.5 < n < 1 and < U < U r . 



and spin systems, providing a general formulation for any 
dimension and any underlying lattice [H-B [ll[ . In this 
paper we have evidenced how the use of this formalism 
leads to the exact solution of the AL-EHM on the Bcthc 
lattice in the presence of an external magnetic field. By 
considering nearest-neighbor repulsion V , there is a tran- 
sition temperature below which a charge ordered phase, 
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FIG. 4: (a) The spin susceptibility as a function of the tem- 
perature for V = 1, n = 0.75, and (a) attractive (b) repulsive 
on-site interaction. 



characterized by a different distribution of the electrons 
in alternating shells, is established for n > 1/z. The on- 
set of the CO phase is signalled by the breaking of trans- 
lational invariance: at T c the values of both the particle 
density and the double occupation become site or, more 
properly, shell dependent Q . The charge ordered phase is 
dramatically affected by the presence of a magnetic field: 
the transition temperature is lessen by a finite magnetic 
field. The CO phase shrinks by increasing h, leading to 
the appearance of large regions where the strength of the 
magnetic field prevents the ordering of the particles. By 
investigating the magnetic properties of the system, we 
found magnetic plateaus at low temperature. Further- 
more, we identified the values of the critical fields h c and 
h s , defining the beginning point of nonzero magnetiza- 
tion and the saturated magnetization field, respectively. 



III. CONCLUDING REMARKS 
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